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SECTION  1 

INTRODUCTION  J 

The  EMP  from  high  altitude  nuclear  explosions  is  traditionally  , 

calculated  in  several  time  phases.  In  the  earliest  time  phase,  which 
covers  times  from  zero  to  about  one  microsecond,  only  time  and  radial 

derivatives  in  Maxwell's  equations  are  significant,  and  only  the  Compton  j 

1 

current  produced  by  previously  unscattered  gamma  rays  is  important.  In 

this  time  phase  the  self-consistent  codes  CHAP  and  HEMP-B  give  quite  accu-  J 

rate  calculations  of  the  EMP.  In  the  time  phase  from  one  microsecond  to 

i 

about  one  millisecond,  scattered  gammas  become  an  important  source  and  the 
transverse  derivatives  in  the  vertical  plane  (of  the  ray  from  burst  to 

I 

observer)  are  not  negligible,  whereas  azimuthal  derivatives  are  still  not 
important.  In  this  case.  Maxwell's  equations  can  be  reduced  to  a set  of 
equations  in  one  time  variable  (the  retarded  time)  and  one  space  variable 
(essentially  the  vertical  coordinate),  in  the  approximation  that  the  gamma- 
ray  front  is  planar  (i.e.,  the  burst  is  assumed  to  be  far  away),  and  the 
curvature  of  the  earth  is  neglected.  This  constitutes  the  far-plane  approxi- 
mation1 upon  which  the  code  LHAP2  is  based.  The  EMP  fields  in  this  time  ' 

phase  are  not  large  enough  to  require  self-consistent  treatment  of  the 
Compton  recoil  electrons,  and  an  approximate  but  fairly  accurate  treatment 
of  the  Compton  current  was  devised  to  increase  the  s^eed  of  the  computations. 

Previously,  the  LHAP  code  has  been  run  only  up  to  retarded  times 
of  100  microseconds.  This  is  approximately  the  time  when  the  EMP  reflected 
from  the  surface  of  the  earth  re-enters  the  source  region,  and  the  code  had 
no  capability  for  handling  the  reflected  pulse.  In  the  present  report  we 


devise  a method  for  including  the  reflected  fields  in  LHAP  so  that  it  can  be 


y.l  ■ 


run  to  later  times.  The  aim  is  to  compute  up  to  about  one  millisecond, 
after  which  other  LI1AP  approximations  break  down. 


The  reflection  of  EM  waves  from  the  earth's  surface  is  an  old 
problem,  and  is  solved  well  in  the  Fourier  domain  by  the  Fresnel  reflection 
factors.3  Several  codes  and  approximate  methods4 * 5have  been  constructed  to 
calculate  the  reflected  EMP.  Most  of  these  methods  Fourier  transform  the 
EMP  time  wave  form,  apply  the  Fresnel  reflection  factor,  and  Fourier  invert 
to  find  the  time  wave  form  of  the  reflected  pulse.  For  the  LHAP  problem 
this  is  not  a convenient  method  for  two  reasons.  First,  the  reflected 
fields  are  needed  at  times  when  the  calculation  of  the  incident  fields  is 
still  in  process;  thus  one  would  have  to  break  the  total  incident  pulse  into 
pieces  of  length  equal  to  the  delay  time  from  source  region  to  ground  and 
back,  and  reflect  each  piece.  Second,  the  processes  of  Fourier  transforming 
and  inverting  are  time  consuming.  It  is  better  to  transform  the  Fresnel 
reflection  factor  to  the  time  domain  in  the  beginning,  for  example,  by 
obtaining  the  reflected  fields  for  an  impulse  incident  pulse.  This  gives 
a kernel  G(t-t'),  and  the  convolution  of  this  kernel  with  an  arbitrary 
incident  pulse  E^(t)  gives  the  reflected  pulse  Er ( t ) , 

t 

Er(t)  = J G(t-f)Ei(f)df  . 

t * =-oo 

This  method  requires  doing  the  convolution  integral  at  each  time  step.  Now, 
if  the  kernel  is  approximated  as  a finite  sum  of  N decaying  exponentials, 
the  convolution  integral  can  be  replaced  by  a set  of  N first  order  dif- 
ferential equations.  The  solution  of  these  equations  needs  only  to  be 
advanced  one  time  step  per  LHAP  time  step,  and  the  computing  time  required 
is  minuscule. 

This  is  the  method  that  we  shall  use.  It  is  equivalent  to 
approximating  the  Fresnel  reflection  factor  in  the  Fourier  domain  as  a sum 
of  poles  l/(w+ibn). 
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SECTION  2 
LHAP  GEOMETRY 


The  physical  situation  is  shown  in  Figure  1.  The  plasma  gamma- 
ray  front  is  shown  at  two  points  in  time  Ct 2>t ^ J . The  ground  plane  and 
the  dashed  planes  limiting  the  desposition  region  are  also  shown.  Each 
star  depicted  is  at  the  same  point  in  space  relative  to  the  nearer  gamma- 
ray  front.  We  assume  that  the  sources  at  star  (1)  and  time  t^,  are  the 
same  as  those  at  star  (2)  and  time  t^.  The  fields  therefore  must  be 
functions  only  of  the  time  after  the  initial  gamma-ray  front  passed  and 
the  altitude.  In  the  code  LHAP2  the  independent  variables  are  defined 
as 


T = ct  - y = retarded  time  in  length  units 
H = z - (tan0)y  = slant  height  above  ground 


where  y is  the  space  coordinate  in  the  direction  n of  propagation  of 
the  gamma-ray  front,  z is  the  coordinate  directed  up  the  wave  front, 

6 is  the  complement  of  the  incidence  angle,  and  ct  is  the  time  coordinate 
in  displacement  units  (c  is  the  velocity  of  light).  These  coordinates  are 
indicated  in  Figure  2.  There  is  no  variation  in  the  x coordinate  which 
is  out  of  the  paper.  The  retarded  time  is  zero  when  the  gamma  front  passes 
a point,  and  x can  be  regarded  either  as  the  distance  of  a point  behind 
the  gamma  front,  or  as  the  time  after  passage  of  the  gamma  front.  In  our 
discussions  we  will  take  the  former  point  of  view  and  think  of  x as  a 
space  variable. 


Changing  to  the  variables  (x,n)  yields  a set  of  first  order 
partial  differential  equations.  These  equations  have  been  differenced  and 
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Figure  1.  Far-plane  approximation  of  LHAP 


using  LHAP  we  can  calculate  their  solutions  by  integrating  forward  in  the 
variable  t.  Near  the  gamma-ray  front  this  gives  the  correct  solution. 

As  indicated  in  Figure  2,  further  to  the  left  of  the  incident  wave  there 
will  be  a region  affected  by  EM  fields  reflected  from  the  earth.  These 
reflected  waves  will  begin  to  have  influence  at  the  point  where  x = s, 
where  s is  the  distance  between  incident  and  reflected  wave  fronts  at 
the  bottom  of  the  deposition  region.  This  distance  is  given  by  the 
following  formula, 

s = 2(sin0)h  . , (2) 

where  is  the  nominal  altitude  of  the  bottom  of  the  deposition  region. 

Before  relating  how  these  waves  are  incorporated  into  LHAP,  we  will  describe 
how  the  reflected  amplitudes  may  be  calculated. 
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SECTION  3 

THE  REFLECTED  WAVE 


During  the  interval  0<t<s  we  have  a complete  solution  for  the 
fields  at  mesh  points  throughout  the  deposition  region.  As  one  approaches 
the  lower  boundary  of  this  region  along  a wave  front  (x  constant),  the 
currents  and  conductivity  become  negligible  and  the  fields  approach  constant 
values.  The  wave  form  along  the  bottom  of  the  deposition  region  is  assumed 
to  be  that  which  reaches  the  earth. 

We  need  to  calculate  the  reflection  of  this  wave.  Scott6  has 
made  measurements  of  the  dielectric  constant  and  conductivity  of  many  soil 
samples,  and  has  found  that  they  are  functions  chiefly  of  the  water  content 
of  the  soil  and  of  the  frequency  of  a monochromatic  wave.  An  RC  network 
model,  in  which  the  resistances  depend  on  water  content,  was  fitted  to 
Scott's  data  by  Longmire  and  Smith.7  Using  this  information  and  the  Fresnel 
formulae  for  reflection  coefficients, 3 it  is  a straightforward  matter  to 
calculate  the  reflected  wave  form. 

A convenient  way  to  calculate  the  reflected  wave  form  is  by  use 
of  the  impulse  response.  Let  the  incident  electric  field  be  the  impulse 
function  6(t).  The  Fourier  transform  of  this  function  is  unity.  Multiplying 
by  the  Fresnel  reflection  factor  p(ui)  (a  complex  function),  we  see  that  the 
impulse  response  is  p(oj)  in  the  frequency  domain.  In  the  time  domain,  the 
impulse  response  is  the  Fourier  inverse, 

00 

G(t)  = ± fp(u>)e~iu)t  dt  . (3) 
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If  the  reflection  factor  is  approximated  as  a sum  of  simple  poles, 

N a 

»<“>  ■ w 

n=l  n 

where  o is  the  infinite  frequency  reflection  coefficient  and  a , b are 
00  n n n 

a set  of  expansion  parameters,  the  impulse  response  becomes 


G(t)  = POT  6(t)  + 0 (t)  £a  e^n1  . 

n=l  n 

where  0(t)  =0  for  t<0  and  0(t)  = 1 otherwise. 


If  we  now  have  an  incident  field  E^(t),  the  reflected  field  is 
the  convolution 


Er(t) 


- I E i 


(t ' ) G(t-t ' )dt ' 


= PooEi^)  + E an  Qn(t)  • (6') 

n=l 

where  we  define 

t 

(^(t)  = J E.(t')  e‘bn(t't,}  dt'  . (7) 

.00 

We  see  that  satisfies  the  simple  differential  equation 

% 

■ Ei«  - bA  • (8) 

Since  the  finite  difference  code  LHAP  advances  the  incident  field  one  step 
at  a time,  it  is  a simple  matter  to  advance  the  at  each  time  step  by 

use  of  Equation  8,  and  then  to  obtain  Er(t)  from  Equation  6' . This  is 
much  easier  than  doing  the  convolution  Equation  6 at  each  time. 


The  incident  wave  at  the  ground  actually  has  two  components  of 
electric  field,  one  in  the  plane  of  incidence  (E^)  and  one  perpendicular  to 
it  (Ex) . The  two  components  have  different  reflection  factors,  each  of  which 
have  to  be  fitted  by  the  form  Equation  4.  The  time  delay  for  the  wave  to 
travel  from  the  bottom  of  the  deposition  region  to  the  ground  and  back  must 
be  taken  into  account. 

The  method  used  here  of  converting  a convolution  integral  into  a 
set  of  differential  equations  has  been  used  previously  by  Longmire  and 
Koppel . 6 


The  method  for  determining  the  expansion  parameters  , an,  b^ 
is  described  in  the  Appendix. 

Figure  3 shows  a typical  incident  electromagnetic  wave  from  LIIAP 
-4 

(up  to  10  second)  and  its  reflection  from  the  earth.  The  reflection  is 
shown  without  time  delay,  i.e.,  for  an  observer  just  above  the  ground.  For 
this  example  we  assumed  a soil  with  ten  percent  water  content  and  an  incidence 
angle  of  30°.  For  the  LHAP  application,  we  need  to  carry  the  reflection  to 
later  times,  of  course. 
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SECTION  4 

MODIFICATION  OF  LHAP 


The  reflected  wave  will  travel  up  into  the  deposition  region  where 
it  will  interact  with  the  ionized  medium  that  lingers  there.  A portion  of 
this  wave  will  be  reflected  by  the  ionized  medium  and  converted  again  into 
a downward- going  wave,  which  will  add  to  the  wave  generated  by  the  late 
Compton  current  sources.  The  sum  of  these  downward-going  waves  will  again 
be  reflected  from  the  ground,  etc.  LHAP  will  calculate  all  of  these  effects 
if  its  boundary  condition  at  the  bottom  of  the  deposition  region  is  modified 
to  include  the  total  fields  reflected  from  the  ground.  That  is,  we  must 
take  the  downward-going  fields  at  the  bottom  of  the  LHAP  mesh  (h  ^ ) , reflect 
them  from  the  ground,  delay  them  in  time,  and  insert  them  as  upward-going 
fields  at  the  bottom  of  the  mesh. 


From  Figure  3 we  see  that  the  earth-reflected  wave  varies  rapidly 
with  time  at  early  times.  In  order  to  handle  these  short  wavelengths  without 
requiring  a fine  n mesh,  we  shift  to  a coordinate  system  fitted  to  the 
reflected  waves.  Figure  4 shows  a reflected  gamma  front  with  propagation 
direction  n'  and  Cartesian  coordinates  y'  and  z';  the  x'  axis  is 
« again  out  of  the  paper.  We  define  the  independent  variables 


» 


x*  = ct  - y'  * retarded  time  for  reflected  wave 
H'  = z’  ♦ (tanb)y'  « slant  height  above  ground 


(9) 


Note  that  t'  is  zero  on  the  reflected  gamma  front,  so  that  there  will  be 
no  reflected  fields  for  t'<0.  The  high  frequency  reflected  fields  arrive 
just  after  t'  * 0,  so  that  a fine  T * mesh  will  be  required  here.  Later 


the  t'  mesh  can  be  expanded,  since  the  reflected  fields  change  more  slowly. 
A fine  n'  mesh  is  not  required.  This  is  the  reason  for  using  the  retarded 
time;  all  of  the  rapid  variations  are  put  into  this  variable  and  removed 
from  the  other. 


Equations  9 differ  from  Equations  1 only  in  the  sign  of  the  tan9 
term.  In  the  earlier  LHAP  writeup2  a parameter  b is  defined  as  tan0.  The 
form  of  Maxwell's  equations  and  the  field  variables  of  that  report  for  the 
incident  wave  frame  can  be  immediately  adapted  for  the  reflected  frame  by 
replacing  b by  -b.  Thus,  in  terms  of  the  Cartesian  field  components 
along  the  x',  y'  and  z'  »xes,  the  modified  LHAP  will  use  the  field  vari- 
ables E , , B , , and 

y”  y' 


H . = E . + B . + bB  , 
x'  x'  z'  y' 


Q,=E,-B,+t-Bi 
Xx ' x ' z ' by' 


H , = E , - B , + bE  , 
z'  z'  x'  y' 


Qz«  = Ez'  + Bx'  " b Ey' 


(10) 


These  equations  replace  Equations  1-14  of  Reference  2.  The  inverse  of  the 
relations  (10)  and  Maxwell's  equations  for  all  of  the  field  quantities  are 
given  by  Equations  1-15  and  1-16  of  Reference  2,  with  b replaced  by  -b. 

We  shall  not  rewrite  these  equations  here,  but  note  some  useful  properties. 
First,  it  is  clear  that  the  new  Cartesian  field  components  are  related  to 
the  old  ones  by  a rotation  through  angle  20  about  the  x (or  x')  axis. 
Thus,  for  the  electric  field  components. 


Ey,  = (cos20)Ey  + (sin26)Ez  (11) 

E , = - (sin20)E  + (cos20)E 
z y z 
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and  similar  relations  hold  for  the  magnetic  field  components.  For  the  H's 
and  Q's,  one  can  then  show,  with  the  help  of  some  trigonometry,  that 


H , = H^  (=0  everywhere) 


Q . = 2E  - Q = H + . B 

x'  xx  x sin20  y 


, = H^  ( = 0 below  source  region) 


Qz,  - 2Bx  - Qz  = - Hz  - ^ 


(12) 


Again  the  Q's  represent  transverse  waves  and  the  H's  longitudinal  waves. 


In  the  regions  outside  the  source  region,  both  of  the  longitudinal 
waves  vanish,  and  the  fields  are  the  superposition  of  transverse  waves 
propagating  downward  (in  the  direction  n)  and  upward  (in  the  direction  n'). 
For  each  of  the  four  transverse  waves  (two  down,  two  up),  the  electric  and 
magnetic  fields  are  equal  but  perpendicular  to  each  other.  We  shall  specify 
the  amplitudes  of  these  four  waves  by  giving  the  x (or  x')  component  of 
the  electric  or  magnetic  field  in  the  down  (d)  or  up  (u)  cases.  In  terms 
of  the  LHAP  field  variables  these  amplitudes  are: 


Exd  = K =-V/si"2e 
Bxd  = 2 = Eyi/sin2e 

(13) 

E = Q . = B /sin29 
xu  2 Xx'  y 

B = i Q . = - E /sin20 
xu  2 xz'  y 

These  amplitudes  will  be  used  in  applying  the  boundary  conditions  at  the  top 
and  bottom  of  the  source  region. 
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In  Figure  2,  the  fields  in  the  region  between  the  incident  front 
and  the  reflected  front  are  to  be  calculated  by  the  previous  version  of 
LHAP.  On  the  reflected  front  t'  =0,  the  coordinate  system  and  the  field 
variables  are  changed  to  the  new  system,  which  is  retained  for  the  duration 
of  the  LHAP  calculation.  Fields  on  the  reflected  front  are  put  in  as  initial 
values  for  the  fields  in  t',  using  Equations  11  and  12. 

We  now  discuss  the  boundary  conditions  to  be  applied  at  the  top 

(h  ) and  the  bottom  f h - ) of  the  n'  mesh.  At  the  top,  the  reason- 
v max  v mm'  r 

able  boundary  condition  is  that  there  are  no  downward-going  waves.  This 
condition  is  met  by  setting  B^,  = E^,  = 0 at  that  boundary,  according  to 
the  first  two  of  Equations  13.  This  parallels  the  procedure  used  in  the 
original  LHAP,  where  = E^  = 0 at  the  bottom  of  the  mesh  to  impose  the 
boundary  condition  of  no  upward-going  waves. 

At  the  bottom  of  the  n'  mesh  (h  . ) , we  pick  out  the  ampli- 

mxn 

tudes  E^  and  B ^ of  the  downward -going  waves,  using  the  first  two  of 
Equations  13.  Each  of  these  amplitudes  is  reflected  from  the  ground,  using 
the  procedure  developed  in  Section  3.  The  reflected  amplitudes,  when  pro- 
perly delayed  in  time  to 


t'  c = t!  + 2(sin0)h  . 
ref  me  v mm 


give  the  values  of  E and  B 

XU  XU 


These  values  determine  boundary  values 


for  Qx,  and  Qz,  from  the  last  two  of  Equations  13. 


These  boundary  values,  including  B^,  = E^,  * 0 at  the  top,  are 
sufficient  to  determine  the  solution  in  the  source  region.  The  two- sweep 
algorithm  used  in  LHAP  is  again  applicable.  Two  trial  values  for  Byl  and 
Eyl  at  the  bottom  are  assumed  and  the  equations  integrated  upward  in  n' 
until  the  top  is  reached,  thus  obtaining  two  trial  solutions.  The  true 
solution  is  then  the  linear  combination  of  these  which  satisfies  the  boundary 
conditions  at  the  top. 
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SECTION  5 
SUMMARY 


A fast  method  has  been  developed  for  computing  the  EMP  waves 
reflected  from  the  ground.  The  LHAP  code  has  been  modified  to  accept  the 
reflected  waves,  which  is  necessary  if  LHAP  is  to  compute  beyond  retarded 
times  of  about  100  microseconds.  A short  test  problem  was  run,  and  the 
code  appears  to  be  working  properly.  Longer  problems  will  be  run  in  the 
coming  year,  when  we  apply  it  to  problems  of  interest  in  EMP  environments. 
We  expect  to  calculate  to  times  of  the  order  of  one  millisecond,  by  which 
time  the  prompt-gamma-driven  EMP  should  be  finished. 
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APPENDIX 

FITTING  OF  THE  IMPULSE  RESPONSE 


* 

I 

I 


We  desire  to  fit  the  impulse  response  function  G(t)  as  a sum 
of  decaying  exponentials.  As  discussed  in  the  text  one  method  of  doing 
this  involves  the  inverse  Fourier  transform  of  the  reflection  coefficients. 
A simpler  procedure  utilizes  the  Laplace  transform  since  the  analysis  will 
then  involve  only  real  variables. 

It  is  useful  to  conceptualize  the  problem  as  an  exponentially 
rising  test  wave  incident  on  the  ground, 

Ei  = Eo  exP(st)  * (Ai) 

The  reflection  of  this  wave  will  also  be  an  exponentially  rising  wave 
with  the  same  rate, 

Er  = pEQ  exp(st)  . (A2) 

The  reflection  coefficient  for  this  wave  has  the  same  form  as  that  for  a 
traditional  Fourier  wave  with  s substituted  for  iw.  For  the  case  with 
the  electric  field  vector  E perpendicular  to  the  plane  of  incidence,  the 
reflection  coefficient  is  given  by 

V2  2 

s(es+47roc)  - s sin  6 

i j 2 2 

s cos0  + \ s(es+4TTac)  - s sin  0 
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When  E is  parallel  to  the  plane  of  incidence. 


(E.S  + 47TOCJ  COs6  - V S (SS+47TQC)  - s^  sin2Q 
(es+47TOc)  cos6  + V s(es+47iac)  - s2  sin20 


In  these,  6 is  the  incidence  angle,  e and  a are  the  dielectric  coef- 
ficient and  conductivity  of  the  earth  respectively,  and  c is  the  speed  of 
light.  Note  that  these  are  both  pure  real  expressions.  The  term  (£s+4ttoc) 
is  the  admittance  of  the  soil  in  cgs  units.  It  is  straightforward  to 
evaluate  this  term  when  the  earth  is  approximated  as  a RC  network,  as  was 
done  by  Longmire  and  Smith.7 


Instead  of  inverse  Laplace  transforming  these  reflection  coef- 
ficients, we  transform  the  response  function  and  make  the  fit  in  the  s- 
domain.  The  Laplace  transform  of  Equation  5 is 

N a 

G(sj  = ow  > Z W 

n=l  n 

We  chose  a p and  a set  of  a and  b so  that  this  function  is  a good 
® n n 6 

approximation  of  the  reflection  coefficient.  p1  and  Pu  must  be  separately 
approximated. 


The  fitting  procedure  is  as  follows.  Let  t . and  t be 
a r min  max 

the  minimum  and  maximum  times  of  interest.  In  the  case  of  LHAP  these  are 


-9 

t . » 10  sec 

min 

t r*  10"3  sec 
max 


(A6) 


Then,  from  the  Fresnel  reflection  factors  for  the  soil  and  angle  of  interest, 

calculate  p(s)  over  the  range  s . £ s £ s , where 

v ' min  max 


MMtttk :A  * T’'  vW; 


IV 


(A7) 

p(sj  will  be  a decreasing  function  of  s,  but  will  not  fall  as  fast  as 

1/s.  Now  choose  a set  of  frequencies  b^  to  cover  the  range  of  s.  Our 

experience  indicates  that  one  frequency  per  decade  will  yield  an  adequately 

accurate  approximation.  Next,  determine  the  and  set  of  a^'s  such 

as  to  give  the  least-squares  best  fit  to  p(s)  using  Equation  A5  at  several 

values  of  s;  say,  s = 103,  3 * 103,  104,  ...,  109.  Finally,  calculate 

the  fit  at  many  values  of  s and  compare  with  p(s).  Better  fits  will  be 

obtained  by  using  more  terms  b and  a . 

6 n n 

In  the  following  table  and  figure  we  present  results  of  a fit  to 
our  example  problem  with  a 30°  incidence  angle  and  a soil  with  a 10  percent 
water  content.  Because  of  the  conductivity,  p(s)  goes  to  unity  for 
s = 0.  Although  we  do  not  need  to  fit  for  s < 103  for  the  LHAP  applica- 
tion, we  note  that  our  procedure  gives 

~ N a 

G(°)  = +Lb^  = °*99968  • (A8) 

n=l  n 

which  is  in  excellent  agreement.  At  the  opposite  extreme,  the  infinite 
frequency  reflection  coefficient  involves  only  the  dielectric  properties 
of  the  soil.  Assuming  the  infinite  frequency  relative  dielectric  coef- 
ficient to  be  5 (see  Reference  7),  p^  has  an  analytic  value  of  0.432, 
somewhat  below  our  fit  value.  Our  procedure  has  lumped  the  responses 

_9 

faster  than  10  sec  with  the  delta  function.  In  the  figure  we  also  show 
the  response  function  with  25  exponentials  as  the  dashed  line.  The  close 
agreement  of  these  curves  makes  it  clear  that  we  have  obtained  an 
excellent  fit. 


» 

* 


Table  1 . 

Expansion  parameters  for  the  reflection 
with  E-field  perpendicular  to  plane  of 
has  a water  content  of  ten  percent,  and 
is  30°. 

coefficient  of  wave 
incidence.  The  soil 
the  incidence  angle 

i 

a^sec"1) 

bi (sec-1 ) 

ai/bi 

1 

2.0257 

103 

0.00203 

2 

2.8295  x io 

104 

0.00283 

3 

1.5070  x 103 

105 

0.01507 

4 

3.5470  x io4 

106 

0.03547 

5 

1.3212  x io6 

107 

0.13212 

6 

1.7462  x 107 

108 

0.17462 

7 

1.0890  x 108 

109 

0.10890 

Delta 

Function  Coefficient 

Poo 

0.52864 

Sum: 


0.99968 


Time  (sec) 


Figure  5 


Response  function.  The  solid  lines  are  the  exponentials  of 
Table  1.  The  dashed  line  is  the  response  function  using 
25  exponentials. 
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